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ABSTRACT 

We show that in 2 + 1 dimensional Quantum Electrodynamics an external magnetic 
field applied to a finite density of massless fermions is screened, due to a 2 + 1- 
dimensional realization of the underlying 2-dimensional axial anomaly of the space 
components of the electric current. This is shown to imply screening of the magnetic 
field, i.e., the Meissner effect. We discuss the physical implications of this result. 
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Quantum Electrodynamics in 2 + 1 dimension has attracted recently consid- 
erable theoretical attention, due to both theoretical and experimental motivation 
from condensed matter systems: on the one hand, several startling features of 
planar quantum mechanics, such as fractional statistics© and C hern- Simons field 
theory© have started revealing themselves; on the other hand both the Quantum 
Hall Effect and High-T c superconductivity are crucially related to the effectively 
two-dimensional nature of the systems where they have been discovered. Indeed, 
the possibility of excitations carrying arbitrary statistics, effectively described by 
coupling of the charged excitations to a Chern-Simons field theory, is at the basis 
of the most succesful theory of the Quantum Hall effect, the Laughlin theory®, 
and is perhaps at the root of high T c superconductivity©. 

The appearance of superconducting behavior in 2 + 1 dimensional systems with 
fractional statistics seems to be a rather robust feature of any theory of fermions 
coupled to Chern-Simons fields®, at least at zero temperature. Because a Chern- 
Simons interaction is always generated dynamically by quantum fluctuations in the 
effective action of 2 + 1-dimensional fermions coupled to external gauge potentials®, 
unless forbidden by symmetry, it is quite conceivable that superconductivity induced 
by the topological structure of the theory which is associated to fractional statistics 
should be a more general feature of 2 + 1-dimensional electrodynamics. 

Here, we shall show that this is indeed true in a well-defined sense: in 2 + 1- 
dimensional electrodynamics with a chemical potential, i.e., with finite fermion 
density, an analogue of the London equation always holds; this (for magnetic fields 
below a critical value) leads to the Meissner effect, i.e., screening of the magnetic 
field, which is a hallmark of superconductivity. 

Rather than analyzing the effective action of fermions coupled to the gauge 
potentials, and discussing possible relations of our result to the usual approaches 
based on an effective Chern-Simons coupling and fractional statistics, we shall es- 
tablish our result by direct computation of the vacuum expectation value of the 
electric current and its curl. We shall show that the the electric current, i.e., the 
space components of the three-dimensional covariant fermionic current of the the- 
ory, satisfies, to first order in the magnetic field, a 2 + 1-dimensional analogue of 
the anomaly equation in two dimensions. This gives to its curl an anomalous, non- 
vanishing vacuum expectation value, proportional to the magnetic field itself, with 
a coefficient that depends on the parity-breaking parameter, i.e., on the chemical 
potential. Thus, quantum fluctuations induce anomalous vacuum currents which 
screen the applied magnetic field, with a screening length depending on the chemical 
potential. We shall establish this result for time-independent gauge fields, as well 
as for time dependent ones at the leading order in the adiabatic approximation, in 
the low temperature limit. We shall also discuss briefly the physical relevance of 
our results. 

We study 2+1 dimensional massless electrodynamics with a chemical potential 
term, defined by the Lagrangian 



tr ^F^F^ + i> (y (z<9 M - - ^ 7 °) V 
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where the coupling e 2 has the dimensions of a mass. We shall work throughout in the 
A = gauge. Our starting point is the observation that the spatial components of 
the fermion current in 2+1-dimensional Minkowski space-time are closely related to 
the axial current of a 2-dimensional Euclidean theory. Indeed, the fermion current 
of the theory (|l|) is 

f = ^yv, ( 2 ) 

where we can represent the 2+1-dimensional Dirac matrices in terms of Pauli 
matrices as Y = io 1 ', (i = 1,2), 7 = a 3 . Consider now the usual Dirac a 
matrices, a 1 = 7°7\ which satisfy a two-dimensional Euclidean Clifford algebra 
{a 1 , a^} = 8^ , and may thus be used to construct a Euclidean Dirac operator 

i p=(id i -A i )a i - 11 

= -HM = -{H + fj). ° } 

The Dirac operator (|3|) is readily recognized to coincide with (minus) the Hamil- 
tonian of the given theory ([!]); notice that the chemical potential term is 

formally identical to a two-dimensional (imaginary) mass term. Recalling that in 
a two-dimensional Euclidean theory 75 = H72)7(2)> ^ follows that the axial current 
associated to the Dirac operator (j^) is 



is (2) = "^7(2)75^ 



(4) 



which is related to the spatial components of the 2 + 1-dimensional Minkowski 
current (@) by 

31(2) = ^°3\ (5) 

where e lJ is the completely antisymmetric tensor in two dimension normalized as 
e 12 = l. 

Now, the current (^) satisfies an axial anomaly equation® 

d tJ l {2) = --B(x) (6) 

where B(x) = eijd l Ai(x). We are thus led to conjecture that for time-independent 
fields, and more generally at the leading order in the adiabatic approximation, the 
spatial part of the current (t, x) at fixed time t should satisfy a similar equation. 
This Eq. would then have the structure of a London equation. A minute's reflection 
reveals that the conjecture cannot be literally true, for the two-dimensional current 
([D and the three-dimensional one (0) do not have the same dimensions, thus the 
space components of the latter cannot satisfy Eq.(H) as is. Also, an equation of 
the form (|||) for the space components of the current signals breaking of parity in 
the 2+1 dimensional theory (just as the presence of a Chern-Simons term in the 
action) and can therefore be induced by radiative correction only in the presence of 
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a parity-breaking parameter. As we shall see, the presence of a chemical potential 
will provide both the parity breaking, and the required dimensionful parameter. 

We shall now proceed to establish our main result, namely, prove the con- 
jecture that the space components of the current (0) satisfy a two-dimensional 
anomaly equation of the form (^). To this purpose, we consider the (^-function 
regularized expression for the vacuum expectation value of the current*© 
Eq.(§: 

(j»(x)) = lim tr [^R(x,x;s)}, (7) 

s— >1 

where R(x, y; s) is the ^-function regularization of the resolvent of the Dirac operator 
which appears in ([!]), i-e., 

i fl 3 lim R(x, y,; s) = 5^ (x — y) 

d „,„A (8) 



s->-l 



where is given by Eq.(|3]). Explicitly we have 

R( X ,v,;s)=^i^fM=^( x \±,\ y) , (9) 

where % p 3 ipk = Xkipk, and the sum, which runs over the full 2+1-dimensional 
spectrum, has been indicated formally as a pointwise trace in the last step, which 
in the x — > y limit we shall henceforth denote with Tr , for short. 

We may now compute the curl of the regulated vacuum expectation value of 
the current (0): 



(e ij diJ j ) = lim Tr 



l 



a j e ij di ' 



(10) 



In order to determine the functional trace in Eq. ([10]) we diagonalize the Dirac opera- 
tor in the adiabatic approximation in terms of the eigenfunctions of the Hamiltonian 
H, by compactifying time over a period T, which we shall let eventually to infinity 
in the zero-temperature limit, and introducing antiperiodic boundary conditions 
over T.(0) The eigenfunctions and eigenvalues of i p 3 (|8|) are then given by 

l p 3 1p nk = Xnk^nk 

(2k + l)n . . 
X n k = y (a n + n) (11) 

tp nk (x y t) = ^e lXnk ^i) n (x) 
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in terms of the so-called Floquet indices©* a n of the instantaneous eigenfunctions 
and eigenvalues of the Hamiltonian (|3]): 



H(x;t)ij) n (x,t) = E n (t)if) n (x,t) 



a, 



1 f T 

- J dtE n {t). 



(12) 



Of course this is the exact spectrum in the case of time-independent gauge fields, 
when a n = E n . 

Substituting Eq.s(|TT|), ([T2|) in the expression ([10|) for the curl of the current we 

get 



(e ij diJ j ) = - lim 



J n i 



tr 101*75 



^ n {x)ijjl{x) 



n k=-oo 

whence, in the particular case of time-independent fields 



(2/c+1)tt . . 

rp "T~ C^ri ~~ r" /-^ 



(13) 



(e ij diJ j ) = lim 



2K tr 



^n75(^)^n(x) 



(2fc + l)7T 



+ E n + fi 



(14) 



= 2 lim Tr j 5 H (D + fx + H) 

s— >1 



where in the last step D = — z J^-, 75 = 03, and we have reverted to the functional 
trace notation which now indicates both summation over the eigenfunctions ( 2fc + 1 ) 7F 
of D and summation and integration over the full spectrum of the Hamiltonian H, 
Eq.@. It is easy to see that Eq.(|l4[) holds in general in the time-dependent case (at 
leading adiabatic order) for the time-average of the curl, provided the eigenvalues 
E n are replaced everywhere by the Floquet indices a n (|i~2|). 

Once the trace over eigenfrequencies of the Dirac operator i p 3 has been re- 
duced to the spectra of the operators D and H according to Eq. ([b|), we may 
proceed to compute the latter. We shall make first use of the fact that both the 
operators H and D have an antisymmetric spectrum, as proven by the existence 
of an operator which anticommutes with each of them, namely 75 and the time- 
reversal operator T, respectively: {75, if} = {T, D} = 0. Using the antisymmetry 
of D first, we get 



(e ij diJ j ) = 2 lim Tr 75 tf (fi + H) (/i + #) ~ D ' 



(15) 



We may now separate the functional traces over D and H by taking a Mellin 



transform, whereby Eq . (|15"1) becomes 

dt 



-2 lim 

S->1 



T(s) 



t 



s-l 



Tr 



D 



-tD 2 



Tr 



H 



K^H + H^e^+V 2 }, (16) 



where Tr d and Tr h denote the traces of the respective operators, which are to 
be evaluated, as usual, by Wick- rotating to Euclidean space and letting T — > 00, in 
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order to enforce the correct (Feynman) boundary conditions. Upon Wick rotation 
the spectrum of H, as well as the chemical potential \i become purely imaginary, 
thereby ensuring convergence of the t integration, and allowing the Mellin transform 
that has lead to Eq.(|l6]), at least if s is kept large enough. 

The trace over D may now be computed explicitly in the large T limit: 



lim Tr d 









2tt J_ 



-tk 2 



1 



(17) 



Using this result in Eq.(|16]) and inverting the Mellin transform we get 

(J'diJ') = - lim Tr ^75 {fiH + H 2 ) + H\}~ s . (18) 

Before computing the trace over the spectrum of the Hamiltonian we may further 
simplify Eq.(|T8"D, using the antisymmetry of the spectrum of H: 



(e ij diJ j ) = lim Tr 7s \H\ \H 2 - a 2 ] 



(19) 



where Tr denotes the trace over all eigenfrequencies of H such that the respective 
eigenvalue A satisfies A > \i (i.e., such that the denominator in Eq.([L9D is positive 
definite). 

The functional trace may now be computed explicitly, observing that in position 
space its diagonal matrix element is given by 



lim (a; I 



Tr<y 5 \H\ [H 2 -^ 2 ] s \x) 



— k^k 2 + <y 5 B(x) \k 2 + ^B(x) - 



■ — \li\B(x) + 



B 

J 2 



(20) 

is as in Eq.(fJ), and as in the standard functional derivation of the 



where B( 

anomaly, ni^P we used the fact that 

H 2 = DiDi + 75 B(x) 



Dj 



-idi + Ai 



(21) 



while, since the trace as s — > is dominated by its large-momentum behavior, the 
spectrum of Di may be approximated by that of 

Using the expression of the trace (|2^) in Eq.(^9|) we get finally 



(e i id i J j (x)) = -^-\ f i\B(x) + 
Ztv 



B 
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(22) 



which is our main result, and proves our conjecture: the vacuum-expectation value 
of the curl of the current is anomalous, and, to first order, proportional to the 
magnetic field. Notice that this a purely quantum effect, in that it pertains a 
vacuum expectation value of the electric current, which is generated by quantum 
fluctuations. 

If we neglect higher order corrections under the assumption that the magnetic 
field is not too strong, Eq . (|22|) is immediately recognized as the London equation, 
which, together with the 2+1-dimensional Maxwell equation 

r = — e ^djB{x) (23) 
leads to the screening equation in the plane 

AB(x) = ^B(x), (24) 
i.e., to the Meissner effect with screening length 




A=W^TTT- (25) 



Let us now discuss the physical implications of this result. If a layered physi- 
cal system may be considered as effectively two-dimensional, and we take the La- 
grangian (|1|) as the effective theory of its charged excitations, with the parameter 
e 2 related to the typical distance between the planes which build up the three- 
dimensional bulk of matter, then a constant magnetic field applied to that material 
is screened, provided in the given region \B\ < \n\ 2 . For example, for high-T c layered 



|2 

superconductors typical values aretLS* e 2 ~ 20 eV. With a chemical potential \i of 
order of a few eV the screening length is A ~ 10 3 A, as in ordinary superconductors. 
When the magnetic flux density exceeds fx 2 , the approximation leading from (|2^ ) 
to ( |2~3| ) breaks down; if we take this as an estimate of the critical flux density we 
get B c ~ (10 3 A)" 2 , i.e., B c ~ 10 3 gauss, as in usual superconductors. 

Even though the present discussion has been given in the low-temperature limit, 
this assumption has been used only in the computation of the trace over Matsubara 
frequencies, Eq.s (|IJ])-([IB|), which has been performed in the limit T — > oo. It ac- 
tually turns out that the first order correction to these results, which is of order 
vanishes, thus suggesting that our results persist up to fairly high temperature (no- 
tice that fi ~ 10 4 °K). Also, even though the chemical potential has played a crucial 
role in our discussion, one may notice that beyond the leading adiabatic approxi- 
mation (to which we have confined ourselves) the spectrum ([II]) also receives a 
Berry phase contribution which enters Eq.(|TT|) in a way formally identical to the 
chemical potential. This suggests that by retaining such contributions our result 
may be generalized to zero density. 

Finally, it is interesting to observe that the induced anomalous currents which 
screen the magnetic field according to Eq. (^2]) are transverse, as those which would 
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be generated by a Chern- Simons term in the effective action of the theory, a char- 
acteristic feature of the Hall effect. This suggests the possibility of applying the 
present results to the Hall system. These and related problems are currently under 
investigation. 

In sum, we have proven that the electric current of fermions in 2 + 1 dimensions 
has an anomalous curl induced by quantum fluctuations, and related to the chiral 
anomaly of an underlying two-dimensional Euclidean theory. This curl obeys the 
Meissner equation and thus it induces screening of the magnetic field. This provides 
a firmer footing to the observation that topological and superconducting behavior 
seem to be rather robust features of gauge theories in 2 + 1 dimensions. It would be 
interesting to try to apply the present considerations to the realistic lattice systems 
which display high-T c superconductivity: perhaps this approach may provide some 
insight in the sought-for novel mechanism which is at the origin of these phenomena. 
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